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We present our progress report on 1+1+1 flavor QCD+QED simulation at the physical point. Cal- 
culations are carried out with 2+1 flavor QCD gauge configurations generated by the PACS-CS 
Collaboration. The dynamical QED effect and the up-down quark mass difference are incorpo- 
rated by the reweighting technique. We also discuss some physics results. 
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1. Introduction 



It has been a long-standing goal to achieve lattice QCD simulations at the physical point. Re- 
cently in [|l|, ^, we have reported the realization of this goal for the 2-1-1 flavor QCD case, where we 
have assumed isospin symmetry and have taken degenerate up and down quark masses, m„ = md. 
However, in nature isospin symmetry is broken by the light quark mass difference, m„ ^ m^, and by 
the difference of the light quark electric charges, 2„ ^ Q^. The isospin symmetry breaking causes 
mass splittings in isospin multiplets of light hadrons, e.g. mj^ — m/f±, m„ — nip. The magnitude of 
splittings are not large and yet important since e.g., it is this difference which guarantees the stabil- 
ity of proton. Because measurements of physical observables of lattice simulations are improving 
beyond 1 % accuracy, we can, and should, take into account the isospin breaking effects. Thus, our 
next target should be l-i-l-i-l QCD-i-QED lattice simulation at the physical point. 

A pioneering QCD-i-QED lattice simulation was carried out in |Q, using quenched approxi- 
mation both for QCD and QED. In this work, the quark fields were explicitly coupled to the photon 
field which was superimposed on gluon fields generated in a conventional Monte Carlo simulation. 
Several studies followed in quenched QED ||5|]-[|^], basically using the above method for including 
electromagnetic effects. The splittings of isospin multiplets of hadrons obtained in these quenched 
QED studies qualitatively matched the experimental ones. Nonetheless, dynamical QED calcula- 
tions were required for convincing results. The possibility of including dynamical QED effects 
were first investigated in [ p^ ] by computing the quark determinant ratio of QCD-i-QED to QCD on 
small lattices with the reweighting technique. Since the stochastic method used for evaluating the 
reweighting factor required many random noises even on a small volume, and scaled as with 
lattice volume W , the calculation on realistically large volumes was expected to be difficult. For 
this reason, most previous attempts were made in QCD-i-quenched QED. 

In this report, we present our preliminary results for l-i-l-i-l full QCD -i- full QED simulation 
at the physical point. The procedure is basically the one developed in [||, [It]]. The 2-1-1 flavor 
QCD configurations on a 32^ x 64 lattice near the physical point generated by the PACS-CS Col- 
laboration are used as the base set for QCD-i-QED. The lattice action is the nonperturbatively 
(^(c?) -improved Wilson quark action and the Iwasaki gauge action at j8 = 1.90. The lattice spacing 
is a O.lfm and the volume is (3fm)^ box. The reweighting technique is used for incorporating 
dynamical QED effects | [T7[ ] and the up-down quark mass difference, and for adjusting the quark 



masses to the physical point P, 18, 19]. 



2. Test calculation of reweighting factor for electric charge 

In order to see the level of difficulty of calculating the reweighting factor for dynamical QED 
effects on QCD configurations with a stochastic method, we investigate the noise dependence of the 
reweighting factor for the strange quark electric charge from eQs = to -^^f^, where ^phys = \j 
and Qs = — |. QCD configurations were generated on a 32^ x 64 lattice [^. Photon fields are 
generated with a non-compact gauge field formulation treated with an appropriate gauge fixing and 
a treatment for zero modes as in [|3|, ^, Employing random noises generated according to the 
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Figure 1: Noise dependence of exponent of re weighting factors for the strange quark electric charge from 



to (ephys2i)/100, where gphys — y/A-n/l'il, Qs — —1/3. Black filled circles show the exponents, 
^ — l)Tji, and have large fluctuation. The even part of the exponents for the electric charge (red 
crosses) is stable. Fluctuation comes mainly from the odd part (blue crosses). 



Gaussian distribution, T]; (/ = 1, . . . ,A^,j), the reweighting factor is evaluated as. 



1 



DetlV= lim _ye-^' (^ -1)''', W 



T-)/£phys6i „ 

^\ 100 ' 

D(0, K,) 



(2.1) 



where D{eQs, Kg) is the Dirac matrix with an electric charge eQs and a hopping parameter K^, W is 
the ratio of the Dirac matrices. 

Figure |] shows the exponent of the reweighting factor riJ{W^^ — 1)t]; for 10 noises labeled 
by / = 1,2, • • • , 10. We see a strong noise dependence (see black filled circles) even though the 
electric charge is turned on at only 1% of the physical value. It is important to observe that the 
main fluctuations come from the ff{e) odd part in the electric charge (blue crosses in Fig. |T]). The 
ff{e'^) even part (red crosses in Fig. |l|) is quite stable. The odd part is not generally forbidden by 
symmetries of the QCD+QED mixed system, so cannot be removed by hand. 

To reduce (or cancel) the {e) odd part, we estimate the reweighting factor for all three light 



quarks simultaneously for each noise |21], 



DetW 



lim -Lye-^H^-'-^)^r 
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D{q'„K,)D{q'^,Kl)D{qi,K,] 
D{qu, K-„) D{qd, Kl) D{qs, K,] 



(2.2) 



where qf = eQf and q'j^ = e'Qf (/ = u,d,s). Since the total charge for up, down and strange quarks 
vanishes, qu +qd + qs = 0, we can expect the cancelation of the odd part, ^(e) = 0, if the quark 
masses are identical. Even though the degeneracy is broken in nature, we expect that ^ (e) terms 
are reduced and so are the fluctuations of the reweighting factor from them. 

To reduce the noise dependence further, we modify the generation of photon fields as follows. 

1) Photon fields are first generated on a 64^ x 128 lattice which is twice finer than the QCD 
lattice of 32^ x 64, by using a non-compact gauge action with an appropriate gauge fixing 



and with a treatment for zero modes as in pL 8, 20] 
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2) To reduce short distance fluctuations of photon fields on QCD lattice, we average the photon 
fields over all independent paths inside the 2^ cell of the QED lattice, which corresponds to 
the unit-cell of QCD lattice. 

3) Finally, compact link fields coupling to quarks with electric charges e^hynQ are defined 
by exponentiating the averaged photon fields A^;U^= e'^p^ysGA^ 

In this procedure, we neglect the running of the electric charge, i.e., take always gphys = \J~^- We 
numerically checked that the averaging procedure did not cause running of the electric charge 
within error bars, by considering electromagnetic mass splittings between charged and neutral 
pseudoscalar mesons in the flavor non-singlet sector and comparing with the case of conventional 
photon fields. We found that the use of the averaged fields reduced the noise dependence of the 
reweighting factor. This procedure also reduced a large additive quark mass shift for Wilson quark 
action due to QED effects to a smaU one. We use these link fields as photons from now on. 



3. Method for 1+1+1 flavor QCD+QED simulation at the physical point 

We carry out 1+1+1 flavor QCD+QED simulation at the physical point as follows, 

1) use 2+1 flavor QCD configurations near the physical point with the nonperturbatively ff{a)- 
improved Wilson quark action and the Iwasaki gauge action at jS = 1 .90 on a 32^ x 64 lattice 
generated by the PACS-CS Collaboration, 

2) generate photon fields as described in the previous section, then superimpose the photon 
fields on the QCD configurations, 

3) reweight into 1+1+1 flavor QCD+QED action with different hoping parameters Ku^d,s and the 
physical electric charge ephys to incorporate dynamical QED effects, 

4) tune three hopping parameters to the physical point values defined by n^,K^,K^ and Q. 
masses as the physical input. 

For the calculation of the reweighting factor, we introduce a square root trick to cut down on the 
calculation time, DetlV = (DetlV)^, and determinant breakups to reduce the noise fluctuation [§, 
18, |l^]. The number of determinant breakups used is 438. 12 noises are used for each determinant 



breakup. We also use the technique of calculating the effects of the up, down and strange electric 



charges simultaneously for each noise as explained in Eq. ( |2.2| ). The inversion algorithm for Dirac 
matrices is the block solver of [ p^ ] which accelerates by a factor of 3 ^ 4 compared with non-block 
solvers. 



4. Preliminary results 

In this section, we present our preliminary results on 1+1+1 flavor QCD+QED simulation at 
the physical point. So far, the number of configurations analyzed is 35 out of 80. 

Fig. ^ shows the configuration dependence of the reweighting factor from 2+1 flavor QCD 
to 1+1+1 flavor QCD+QED at the physical point. Except for a few configurations, the fluctuation 
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Figure 2: Configuration dependence of the reweighting factor from 2+1 flavor QCD to 1+1+1 flavor 
QCD+QED. 
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Figure 3: Noise dependence of m„-+,m[^+ and m^o normalized by m^- (black filled circles). Red lines 
indicate the experimental values. The horizontal lines indicate the number of noises used in each determinant 
breakup for the reweighting factor. 



is contained within iff (10) due to the introduction of the techniques as described in the previous 
sections. 

For hadron measurements, four source points were used. We tuned the hopping parameters to 
the physical point values using 71+,^'+,^*' and masses as the physical input. Figure || shows the 
hadron mass ratios, and (black filled circles), as a function of the number of noises 

used for each determinant breakup. The ratios reach a plateau around the experimental values (red 
lines) already from small number of noises. This shows that twelve noises for each determinant 
breakup is sufficient for the reweighting factor calculation. Figure |3| also shows that the tuning to 
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Figure 4: Ratio of to propagators which shows the mass difference m^o — + . Black curve is the 
lattice result which is roughly consistent with the experimental values (red curve). Blue curves are fit ones 
with 1 a error band. 



the physical point has been achieved up to error bars. 

Figure ^ shows the ratio of to propagators which is expected to behave as a function of 
time t as, 



{K\t)K^\Q)) 
{K+it)K+{0)) 



(4.1) 



where we assume (m^o —mic+)t <^ 1. The slope indicates the mass difference. The lattice data 
(black curve) and the fit result with 1 a error band 3 .2 1 (57) MeV (blue curve) are roughly consistent 



with the experimental value 3.937(28) MeV [23] (red curve). 

Finally in Table 1, we summarize the up, down and strange quark masses renormalized at /i = 2 
GeV in the continuum MS scheme. The renormalization factor is nonperturbatively determined in 
the Schrodinger functional scheme [^]. Here, we neglect QED corrections to the renormalization 
factor. 

The results presented here were obtained from 35 configurations out of 80 ones and are still 
preliminary. In order to conclude 1+1+1 flavor QCD+QED simulation at the physical point, we 
are now increasing the statistics and also source points for measurements to reduce the statistical 
errors. 



mu 
md 

rris 

(m„ + md)/2 
mu/md 

2ms/{mu + md) 



1.97(67) [MeV] 
4.31(83) [MeV] 
90.32(67) [MeV] 
3.14(72) [MeV] 
0.457(93) 
28.8(6.6) 



Table 1: Preliminary results for renormalized quark masses in 1+1+1 flavor QCD+QED at j3 = 1.90 on 
32^ X 64 lattice. = 2GeV. 
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